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This  report  derives  a  mathematical  model  that  describes  the  minimum- 
weight  design  of  a  horizontal  beam  under  an  invariable,  concentrated, 
vertical  load  moving  slowly  from  one  end  to  the  other.  Beams  of  this 
type,  although  possibly  impractical  in  earthbound  structures,  may 
someday  become  interesting  for  applications  in  a  weightless  environment 
such  as  that  of  a  space  station. 
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PREFACE 

Plastic  design  of  structures  has  received  considerable 
attention  in  recent  years.  To  the  authors'  knowledge,  the 
present  Project  RAND  memorandum  is  an  innovation  in  this 
field  and  also  in  the  domain  of  linear  programming  —  in  the 
former  because  the  beam  design  Is  to  withstand  a  one— parameter 
family  of  load  condltj  ons  rather  than  a  single  load  distribu¬ 
tion,  and  in  the  latter  because  the  program  involves  a  continuum 
rather  than  a  finite  number  of  variables. 

The  basic  mathematical  treatment  of  the  continuous  version 
of  the  problems  presented  in  this  paper  is  contained  in  RM— 2993  PR, 
"A  Linear  Program  of  Prager's:  Notes  on  Linear  Programming  and 
Extensions  —  Part  60,"  now  in  the  process  of  publication.  In 
the  present  paper,  which  is  to  be  submitted  for  presentation  at 
the  Fourth  U.  S.  National  Congress  of  Applied  Mechanics  at 
Berkeley,  California,  June,  1962,  the  mathematical  model  is  derived 
from  physical  considerations  and  solved  in  a  form  readily  accessi¬ 
ble  to  the  structural  engineer. 

It  might  be  added  that  beams  of  the  type  considered  in  this 
paper  may  be  impractical  (i.e.,  expensive  to  manufacture)  for 
earth-bound  structures.  They  are,  however,  important  for  theo¬ 
retical  reasons  since  they  represent,  in  a  sense,  the  ultimate 
in  weight  design.  It  Is  of  course  conceivable  that  such  designs 
might  become  more  important  in  the  future  as  actual  components 
of  nonearth— bound  structures  such  as  space  stations. 
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SUMMARY 

X 

The  main  problem  considered  and  solved  Jm«6M«e9ft£er 
Is  the  minimum-weight  design  of  a  horizontal  I— beam  of  con¬ 
stant  web  and  variable  flange  thickness.  The  beam  is  simply 
supported  at  one  end  and  built  in  at  the  other  and  is  designed 
to  withstand  (in  plastic  flow)  a  concentrated  vertical  load 
of  fixed  intensity  moving  slowly  from  one  end  to  the  other. 

To  introduce  the  raarflrar— frit  the  general  techniques 
employed,  a  few  discrete  versions  of  the  problem  are  first 
presented  and  solved  as  linear— programming  problems. 

j 
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MINIMUM-WEIGHT  DESIGN  FOR  MOVING  LOADS 
1. _ INTRODUCTION 

Plastic  design  of  continuous  beams  and  frames  for  mini¬ 
mum  weight  has  been  extensively  studied  in  recent  years. 

Almost  without  exception,  earlier  studies  dealt  with  structure 
of  given  over— all  dimensions  (widths  of  bays  and  heights  of 
stories)  that  are  subjected  to  a  single  system  of  concentrated 
loads  judged  to  be  the  most  dangerous  combination  of  the 
groups  of  loads  that  the  structure  has  to  carry.  For  practi¬ 
cal  reasons,  the  beams  and  columns  of  steel  frames  are 
usually  prismatic,  so  that  the  designer  disposes  of  only  a 
finite  number  of  design  parameters  (the  yield  moments  of  the 
beams  and  columns).  If,  within  the  range  of  practical 
structural  sections,  the  unit  weight  of  a  section  can  be 
approximated  by  a  linear  function  of  its  yield  moment, 
designing  for  minimum  weight  can  then  be  formulated  as  a 
problem  in  linear  programming  (see,  for  Instance,  [l]*,p.  84 
ff . ) . 

While  designs  with  beams  or  columns  of  continuously 
varying  cross  section  would  not,  in  general,  be  practical, 
they  are  worth  Investigating  because  they  furnish  the 
theoretical  minimum  of  the  structural  weight,  against  which 
practical  designs  may  be  checked. 

*Numbers  In  brackets  Indicate  items  in  the  list  of 
references  at  the  end  of  the  paper. 
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The  minimum-weight  design  of  a  beam  with  continuously 
varying  cross  section  subjected  to  a  single  system  of  loadB 
has  been  discussed  by  Heyman  [_  2  Consider,  for  instance, 
a  sandwich  beam  with  a  core  of  given  constant  height  h  and 
breadth  b,  and  flanges  of  constant  breadth  b  and  variable 
thickness  t.  Since  the  designer  controls  only  the  thickness 
t,  the  variable  part  of  the  structural  weight  is  proportional 
to  the  integral  of  t  extended  along  the  beam.  Using  a  theorem 
of  Drucker  and  Shield  C  3  J  »  Heyman  has  shown  that  the  abso¬ 
lute  minimum-weight  design  of  such  a  beam  can  be  based  on  a 
deflected  shape  for  which  the  curvature  has  constant  abso¬ 
lute  value.  Figure  1  shows  deflections  of  this  kind  for  a 

beam  that  is  simply  supported  at  one  end  (x  =  0)  and  built 

in  at  the  other  end  (x  =  i).  The  curvature  is  found  to 

change  sign  at  x  =  l/'i/’ZT .  Accordingly,  the  bending  moment 
must  vanish  at  x  =  Since  the  indeterminacy  of  the 

beam  is  removed  by  this  condition,  the  bending  moments  can  be 
obtained  from  the  design  loads  without  reference  to  the  cross- 
sectional  dimensions  of  the  beam.  The  flange  thickness  at 
each  cross  section  is  then  determined  from  the  bending 
moment  at  this  section  in  such  a  way  that  the  yield  limit  of 
the  flange  material  is  reached  at  all  sections  when  the  beam 
is  subjected  to  the  design  loads. 

As  this  example  shows,  an  absolute  minimum-weight  design 
for  a  single  system  of  loads  is  essentially  a  statically 
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determinate  design.  This  feature,  which  greatly  facilitates 
the  design,  is  likely  to  be  absent  when  a  structure  has  to 
be  designed  for  more  than  one  system  of  loads.  To  illustrate 
the  arguments  that  arise  in  the  solution  of  problems  of  this 
kind,  the  sandwich  beam  in  Fig.  1  will  be  designed  for  a 
moving  load  of  constant  intensity  P  that  may  act  at  any 
section  of  the  beam.  Calling  for  the  determination  of  the 
flange  thickness  at  each  station  x,  this  problem  involves 
an  infinity  of  design  parameters  or,  what  amounts  to  the 
same,  a  design  function  t(x).  Since  there  are  as  yet  no 
general  methods  of  solving  continuous  problems  of  this  kind, 
a  discrete  analog  will  be  discussed  in  Sec.  2  as  a  guide 
to  the  solution  of  the  continuous  problem,  which  will  be 
given  in  Sec.  3. 

2.  DISCRETE  POSITIONS  OF  THE  LOAD 

The  general  features  of  the  desired  design  can  be 
explored  in  the  following  way.  Dividing  the  span  £  of  the 
beam  into  n  equal  segments,  suppose  that  only  the  partition 
points  are  eligible  as  points  of  application  of  the  load,  and 
that  the  flange  thickness  t  is  a  segmentwise  linear  function 
of  the  distance  x  measured  along  the  beam.  This  kind  of 
design  is  specified  by  the  n  values  of  t  at  the  points  x  =  i/n 
2£/n,  ...,  (n— l)£/n,  £  ,  the  flange  thickness  being  zero  at 
the  simply  supported  end  x  =  0.  It  is  to  be  expected  that 
the  desired  design  with  continuously  varying  flange  thickness 
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results  as  n  — >  co  . 

If  h  denotes  the  constant  height  and  b  the  constant 
breadth  of  the  core,  and  sq  the  tensile  or  compressive  yield 
stress  of  the  flange,  and  if  the  flange  thickness  t  is 
treated  as  small  in  comparison  to  h,  the  yield  moment  is 

(2.1)  Y  =  so  bht. 

For  a  unit  length  of  the  beam,  the  variable  part  of  the 
structural  weight,  that  is,  the  weight  of  the  flanges,  is 
proportional  to  the  flange  thickness  t  and  hence  to  the 
yield  moment  (2.1). 

For  n  =  2,  a  design  of  the  considered  type  is  specified 
by  the  yield  moments  Y-^  and  Y2  at  the  center  of  the  beam 
and  the  built-in  end.  Only  one  position  of  the  load  and 
one  yield  mechanism  need  to  be  investigated.  With  the 
deflections  and  hinge  rotations  shewn  in  Fig.  2,  the  condi¬ 
tion  that  the  energy  dissipated  in  the  plastic  hinges  equals 
the  work  of  the  load  takes  the  form 

(2.2)  2Y1  +  Y2  =  Pi/2. 

Since  the  yield  moment  vanishes  at  the  simply  supported  end, 
the  average  values  of  the  limiting  moment  in  the  two  halves 
of  the  span  are  Y-j/2  and  (Y1  +  Y2)/2,  so  that  the  variable 
part  of  the  structural  weight  is  proportional  to 


(2.3) 


W  =  (2Y1  +  Y  2)l/h 
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The  yield  moments  Y^  and  Yg  of  a  minimum-weight  design  must 
satisfy  (2.2).  Substitution  of  this  equation  into  (2.3) 
shows  the  variable  part  of  the  structural  weight  to  be 
proportional  to 

(2.4)  W  =  ?£%. 

Thus  there  exist  infinitely  many  minimum-weight  designs 
for  n  =  2. 

For  a  greater  value  of  n,  a  yield  mechanism  is  defined 
by  the  locations  of  its  two  plastic  hinges,  say  x  -  h£/n 
and  x  =  k£/n  with  k  >  h.  From  the  deflections  and  hinge 
rotations  Indicated  in  Fig.  3,  the  condition  that  the  energy 
dissipation  in  the  plastic  hinges  must  not  exceed  the  work 
of  the  load  becomes 

(2.5)  kYh  +  hYk  ^  Ph (k— h) £/n . 

The  variable  part  of  the  structural  weight  is  proportional 
to 

(2.6)  W  =  (Y1  +  Y2  +  ...  +  Yn_x  +  0.5  Yn)i/n. 

For  n  =  4,  for  instance,  the  (nonnegative)  yield  moments 
of  the  minimum-weight  design  must  minimize  the  linear  form 


(2.7) 


W  *=  (Yx  +  Y2  +  Y  +  0.5  Y4)£/\ 


FIG.  5 
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while  satisfying  the  inequalities 


2Yl  +  *2 

2. 

3*1 

£  Pi/2 

4*1 

+  y4 

£  3P lA 

3Y 2  +  2Y3 

^  Pi/2 

4*2 

+  2Y4 

>  Pi, 

4Y5 

+  3Y4 

£  3PiA 

which  are  obtained  from  (2.5)  by  letting  h  and  k  vary  from 
1  to  4  subject  to  the  condition  that  k  >  h. 

In  general,  the  n  yield  moments  of  the  minimum-weight 
design  must  minimize  a  linear  form,  which  represents  the 
variable  part  of  the  structural  v/eight,  while  satisfying 
n(n— l)/2  inequality  constraints.  Expressed  in  structural 
terms,  the  essential  features  of  the  solution  of  this 
linear— programming  problem  are  as  follows  (see,  for  instance, 
14]).  The  yield  moments  of  the  minimum-weight  design  ful¬ 
fill  n  of  the  n(n— l)/2  constraints  as  equations  rather  than 
inequalities.  Each  of  these  equations  corresponds  to  a 
yield  mechanism  that  becomes  critical  for  some  position  of 
the  moving  load;  it  states  that  the  energy  dissipated  In 
the  plastic  hinges  of  this  mechanism  must  equal  the  work  of 
the  load  if  the  mechanism  is  to  be  critical.  The  n  equations 
can  be  linearly  combined  with  nonnegative  factors  in  such 
a  manner  that  the  coefficient  of  each  yield  moment  in  the 
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linear  combination  equals  the  coefficient  of  this  yield 
moment  in  the  expression  for  the  variable  part  of  the 
structural  weight. 

For  n  =  4,  the  use  of  the  equality  sign  in  the  second, 
fourth,  fifth,  and  sixth  relations  (2.8)  furnishes  the  yield 
moments 

(2.9)  Y1  =  5Pi/32,  Y2  =  7Fi/48,  Y3  =  Pi/52,  Yj,  =  5Pi/24, 

which  are  readily  seen  to  satisfy  the  remaining  two  relations 
(2.8).  Moreover,  when  the  four  equations  corresponding  to 
critical  yield  mechanisms  are  combined  with  the  factors  i/12, 
i/24,  i/52,  i/48,  the  coefficient  of  each  yield  moment  in 
the  linear  combination  is  found  to  equal  the  coefficient  of 
this  yield  moment  in  (2.7).  The  yield  moments  (2.9)  there¬ 
fore  represent  the  minimum— weight  design. 

In  Fig.  4,  the  line  OABCD  represents  the  variation  of 
the  yield  moment  along  the  beam;  and  OA'B'C'D'  is  the  image 
of  this  line,  by  symmetry  with  respect  to  the  x-exls.  For 
the  considered  positions  of  the  load  at  x  =  1/ 4,  x  =  i/2, 
and  x  =  5  i/4,  the  bending  moments  are  represented  by  "the 
lines  OAC'EjOBC'D’ ,  and  OCD' ,  respectively.  For  none  of 
these  load  positions  does  the  absolute  value  of  the  bending 
moment  at  any  section  of  the  beam  exceed  the  yield  moment 
at  this  section,  while  for  each  load  position  the  absolute 
value  of  the  bending  moment  attains  the  yield  moment  in  at 
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least  two  sections. 

With  the  yield  moments  (2.9)>  (2.7)  furnishes 
(2.10)  W  =  7P l2 /6k. 

3.  CONTINUOUSLY  VARYING  POSITION  OF  THE  LOAD 

The  solution  of  the  discrete  problem  in  Sec.  2  suggests 
that  for  continuously  varying  position  of  the  load  the  yield 
moments  of  the  minimum  weight  design  are  represented  by  arcs 
such  as  OABC  and  CD  in  Fig.  5*  Denoting  the  abscissa  of  C 
by  x.,  suppose  that  the  load  acts  at  the  section  x  =  xq  <  x1 
and  that  a  plastic  hinge  forms  at  this  section.  The  bending 
moment  at  this  section  then  equals  the  yield  moment  Y(xq), 
and  the  reaction  at  the  simply  supported  end  equals 

(5-1)  R(x0)  =  Y(x0)/xo' 

To  the  left  of  the  load,  the  bending  moments  are  given  by 

(5.2)  M(x;xQ)  =  x  R(xo)  =  x  Y(xo)/xo  ,  0  <  x  ^  xq 

(line  OB  in  Fig.  5).  If  the  arc  CABC  is  concave  with  respect 
to  the  x-axis,  these  bending  moments  nowhere  exceed  the 
yield  moment.  To  the  right  of  the  load,  the  bending  moments 
are  given  by 

(5-5)  M(x;xq)  =  x  R(xo)  -  (x  -  xo)P  ,  xQ  <  x  ^  £ 

(line  BE  in  Fig.  5).  Since  there  must  be  a  yield  mechanism 


RM— 2887 

11 


PIG.  5 


RM-2887 

12 


corresponding  to  the  considered  position  of  the  load,  the 
absolute  value  of  the  bending  moment  must  attain  the  yield 
moment  at  some  abscissa  x  =  x2  ^  x^  without  exceeding  the 
yield  moment  at  any  other  abscissa.  In  other  words,  the 
line  BE  must  be  tangent  to  the  image  CD1  of  the  arc  CD  with 
respect  to  the  x-axis.  Thus,  one  obtains 


(3.4)  Y(x)  +  M(x;xq)  ^  0  ,  xx  <  x  ^  i, 

with  equality  for  one  value  of  x.  With  the  use  of  (3.1) 
and  (3-3)»one  finds 


(3.5)  Y(x)  1  (x  -  xQ)P  -  x  Y(xo)/xo  ,  x±  £  x  <  f. 


To  find  the  abscissa  of  the  point  of  contact  E  of  the  line 
BE  and  the  arc  C'D',  equate  the  derivative  of  the  right-hand 
side  of  (3.5)  to  zero  and  replace  x  by  x^.  Thus,  x2  Is  given 
by 


(3-6) 


F 


0  <  xo  <  X1  » 


where  the  prime  denotes  differentiation  with  respect  to  xQ. 
V/ith  this  value  of  x2,  the  yield  moment  Y(x2)  is  found  by 
setting  x  =  x2  and  using  the  equality  sign  in  (3*5).  As 
xq  varies  from  0  to  x^,  the  point  of  contact  E  moves  from 
C'  to  D'  as  indicated  in  Fig.  5.  Accordingly,  the  arc  C'D' 
Is  uniquely  determined  by  the  arc  OABC. 

The  variable  part  of  the  structural  weight  is  propor— 
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tional  to 

x  0/ 

(3.7)  w  =  J  1  Y(x)  dx  +  f  Y(x)  dx. 

J  o  J 

Let  Y(xo)  be  an  arbitrarily  defined  positive  concave  function 
in  0  £  xq  £  x^  that  vanishes  for  xQ  =  0,  and  x  =  f(xo)  an  arbi¬ 
trary  increasing  function  that  maps  the  interval  (0,x1)  onto 
(x1#£).  When  the  second  integral  in  ( 3 - 7 )  is  transformed  by 
the  use  of  the  inequality  (3*5)  and  the  change  of  variable 
x  =  f(xo),  the  relation 

(3.8)  w  ^  r*1  C*(*o>  t1  -  +  (f  -  xo)f'P]  dxo 

1  o  o 

is  obtained.  To  render  the  lower  bound  (3*6)  for  the  struc¬ 
tural  weight  independent  of  the  choice  of  Y(xo),  set 

(3-9)  ff  =  xc, 

or 

(5-10)  f  +  c2  , 

where  c  is  a  constant  of  integration. 

Since  f  is  to  map  the  interval  (0,x^)  onto  (x^,i!),  f  satisfies 
(3.11)  f(0)  =  x2  ,  f(x1)  =  l  . 

It  follows  from  (3. 10)  and  (3. 11)  that 


(3.12) 


c  =  xx  =  1//T  . 


RM-2887 

14 


and 

(5.15) 


In  view  of  (5-9)  and  (3. 13),  the  inequality  (3.8)  reduces  to 
(3.14)  W  'l  J(1  -V?  +  log(l  +</?))P£2  =  0.1168P£2. 


To  obtain  the  function  Y(xq)  for  which  W  equals  this 
lower  bound.  Identify  the  function  (3-14)  with  the  right-hand 
side  of  (3.6)  to  obtain 


(3.16)  Y(xo) 


Pxo  {k  “  log  t(xo 


where  k  Is  a  constant  of  integration.  To  determine  this  constant, 
note  that  Y(xx)  evaluated  from  (3-18)  must  equal  Y(x1)  obtained 
by  substituting  xQ  =  0  and  x  =  x^  into  (3«5)  and  using  the 
equality  sign  in  this  relation.  Thus,  the  value  of  k  is 

k  =  ^(1  +  lo6  1  g1^)  • 

With  f  =  x2,  it  follows  from  (3.13)  that 

(3-18)  xo  =1/x2  -  x2. 


and  from  (3-5)  and  (3«l6)  that 
(3.19) 

r  /“5 - 2  /"S - ? 

Y(x2)  =  P  |x2  -i/x2  ~  xt  ~  x2^x2  ”  X1  k 

-  log[(x2  +  i/x2  ~  }»  xi  <  x2  ^ 
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where  k  Is  given  by  (3-17) •  The  arcs  OABC  and  CD  in  Fig.  5 
represent  the  functions  (3.16)  and  (3*19)* 

It  is  interesting  to  note  that  the  abscissa  defined 
by  (3-12)  is  identical  with  the  abscissa  at  which  the  curvature 
in  Fig.  1  changes  sign.  An  analogous  identity  was  found  to 
exist  for  the  doubly  built-in  beam,  which  was  treated  along 
similar  lines,  but  no  general  reason  for  this  has  been  dis¬ 


covered  . 
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